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We study the electroweak phase transition by lattice simulations of an effective 3-dimensional theory, for a 
Higgs mass of about 70 GeV. Exploiting a variant of the equal weight criterion of phase equilibrium, we obtain 
transition temperature, latent heat and surface tension and compare with Mh ~ 35 GeV. For the symmetric 
phase, bound state masses and the static force are determined and compared with results for pure SU (2) theory. 



1. Introduction 

There is an old belief that the electroweak stan- 
dard theory possesses a first order phase transi- 
tion at some temperature of the order of 
the W mass. This phase transition has become 
subject of intensive studies in the last years, in 
particular its dependence on the mass of the so 
far elusive Higgs boson. One motivation was the 
phenomenological interest in baryon asymmetry 
generation at the electroweak scale. The tran- 
sition has to be strong enough, both in order 
to accomplish a sufficient rate of baryon gener- 
ation during the transition and to prevent the 
wash-out of baryon number after it is completed. 
The present quantitative understanding of pos- 
sible mechanisms as well as experimental lower 
bounds for the Higgs mass make this unlikely 
within the minimal standard model. 

A second reason was the wish to control the 
behaviour of perturbative calculations of the ef- 
fective action. This quantity is the appropriate 
tool of (non-lattice) thermal quantum field the- 
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ory for dealing with symmetry breaking. Infrared 
problems prevent a perturbative evaluation of the 
free energy in the symmetric phase to higher loop 
order. However, the true, non-perturbative na- 
ture of the symmetric phase will be characterized 
by massive W- and Higgs bound states instead 
of massless W gauge bosons. A selfconsistent ap- 
proach to provide masses across the transition, 
e.g. by gap equations, can improve the ability 
to calculate perturbatively the symmetric phase. 
Gauge field condensates are another property of 
the symmetric phase, expected to lower its free 
energy density. 

An independent method is needed to character- 
ize the electroweak phase transition even within 
a pure SU (2) gauge-Higgs version of the theory. 
This model has become a testfield to control the 
validity of perturbative predictions over a broad 
range of Higgs masses. At present, one is inter- 
ested to see whether the first order transition ends 
somewhere around a Higgs mass Mjj 100 GeV. 
Lattice simulations @-|]l^ are not only able to 
describe both phases starting from first princi- 
ples but, moreover, make it possible to put both 
phases into coexistence near the phase equilib- 
rium. Thus one is able to measure directly quanti- 
ties like latent heat, surface tension, condensates 
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etc. quantifying the strength of the transition. 

One approach to lattice calculations of the elec- 
troweak transition is based on an effective 3- 
dimensional Higgs model. It is attractive phe- 
nomenologically because it circumvents the prob- 
lem of putting chiral fcrmions on the lattice. 
Due to dimensional reduction, fermions as well 
as non-static bosonic modes contribute to the ef- 
fective action. In contrast to QCD, dimensional 
reduction should work for the electroweak the- 
ory around and above the transition tempera- 
ture because is small. For the electroweak 
phase transition this approach has been pioneered 
by Farakos et al. (see e.g. [§,|ll). This pro- 
gram aims at exploring the accuracy of dimen- 
sional reduction at various Higgs masses by com- 
paring various parameters of the transition with 
those of 4-dimensional lattice and perturbative 
approaches. Perturbation theory is necessary to 
relate the 4-dimensional continuum theory to the 
parameters of the dimcnsionally reduced theory 
and, finally, to the bare coupling parameters of 
the lattice action. Dimcnsionally reduced ver- 
sions retain the remnant of the temporal gauge 
field (as an adjoint Higgs field) or not (as in 
this work). 

Recently [|j we have presented results obtained 
with the 3-dimcnsional lattice model on the phase 
transition for « 35 GeV. Here we present 
some results of numerical work on the more re- 
alistic Mh ~ 70 GeV and compare them with 
the case of smaller Higgs mass. As expected the 
first order nature has become weaker but is still 
evident. In this talk we put our main empha- 
sis on ways to characterize the phase equilibrium 
at finite lattice size in order to obtain the infi- 
nite volume limit of the transition parameters, 
and on non-perturbative features of the symmet- 
ric phase. The extrapolation to the continuum 
limit will be dealt with in a forthcoming publica- 
tion [|2|. 



2. The Model 

We study the S'[/(2)-Higgs system with one 
complex Higgs doublet of variable modulus. The 
gauge field is represented by the unitary 2x2 link 
matrices Ux,^^ and the Higgs fields are written as 
$x = PxVx [pi = \Tr{^+<^x) is the Higgs mod- 
ulus squared, T4 an element of the group S'[/(2)). 
The lattice action is 

S = /3G^(l-irrt/p)- 
p 

I 

J2{pl + Pnipl-m (1) 

X 

with /3g = 4/{ag'^). In three dimensions the lat- 
tice Higgs self-coupling is /3fl = {h/gl) (Ph/Pg), 
and A3 denote the 3-d continuum gauge and 
Higgs self couplings which are 3-d renormaliza- 
tion group invariants. They are related to the cor- 
responding four dimensional couplings via g^ = 
Tig'' + 0{g^)) and A3 = T(A + 0{g^)). 

String operators like 
E{1) = <i>tU^^f,Ux+f,,f.-Ux+(i-i)ii,f,^x+ifi (2) 

(of extension I) are used to form Higgs and W- 

operators. Actually the 3 d masses (inverse 

correlation lengths) have been obtained from the 
connected correlators between separated "time 
slice" (in 2 + 1 dimensions !) sums of "spatial" 
string operators, which project out the proper 
SU{2) and spin content. The overlap with the 
lowest mass states is improved choosing an ex- 
tension Z = 4 for correlators in the W and also 
partly in the Higgs channel (there also is 
used). A static force is (formally) defined using 
Wilson loops W{R, T) of asymmetric extensions 
2 < i? < r < L/2 (L^ is the lattice size). 

The bare lattice coupling parameters (3h and 
I3r can be translated into a physical tempera- 
ture by a formula which generalizes the tree level 
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relation between the Higgs mass and the bare 
lattice parameters including perturbative correc- 
tions (too complicated to be given here) . We con- 
sider here the quartic coupling Aa/gl = 0.0957. 
According to the tree level based relation A3 / « 
Xjg^ « M^/(8M^) this would correspond to the 
case of Mh = 70 GeV. Corrections to this approx- 
imate relation depend on details of the dimen- 
sional reduction (loop order, the adjoint Higgs 
field integrated out or not). Taking these cor- 
rections into account (4-dimensional renormal- 
ization effects are neglected) this coupling ratio 
corresponds to Mh = 71.8 GeV for = A^j^. 
For comparison, the lower Higgs mass had been 
simulated with g\ = 0.0239, corresponding to 
Mh = 38.3 GeV. 

The lattice Higgs self coupling (3j^, is not fixed 
but runs with Ph- This is important for the mul- 
tihistogram technique to be used. In general, bin- 
ning has to be performed in the two relevant parts 
of the action (corresponding to (3h and Pr) and 
some other observable. Our data come from typi- 
cally 100000 configurations per set of couplings 
(separated by one heat bath step mixed with 
eight Higgs field reflections). Correlation mea- 
surements were separated by 10 iterations. 



3. Localization of the phase transition 

In order to characterize the order of the phase 
transition at Higgs mass AIh ~ 70 GeV and to 
obtain the infinite volume critical hopping pa- 
rameter (3^^ we have simulated lattice sizes 30'^, 
48^^ and 64'^. To control the approach to the 
continuum limit we have chosen gauge couplings 
Pg = 12 and Pa = 16 (not presented here). Lat- 
tice observables to monitor the phase transition 
are the Higgs modulus squared and the link 
contribution to the action -E(l), but other pure 
gauge field quantities like the average plaquette 
or the Polyakov line give 2-state signals, too. 

We have used the Ferrenberg-Swendsen analy- 
sis and present here only results obtained for his- 
tograms of p^. The corresponding Binder cumu- 
lants B^2 are shown for the lattice sizes L = 30, 



48 and 64 in fig. ^ The minima give the respec- 
tive pseudocritical Ph^. From a linear extrap- 
olation with we obtain the infinite volume 
limit Phc = 0.3435433. For comparison, the max- 
ima of the p^ susceptibility give Phc = 0.3435430. 
In the case of a SU{2) Higgs model without t- 
quark this value corresponds to Tc = 156 GeV. 
For comparison, we obtain Tc = 98.7 GeV for 
Mh « 35 GeV. Finiteness and shrinking of the 
Binder cumulant with increasing volume are evi- 
dence for the first order nature of the transition 
at the larger Higgs mass, too. The dip of the 
Binder cumulant reaches 0.6623 in the V —^ 00 
limit clearly different from 2/3. 
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Figure 1. Binder cumulants B^2 for L = 30, 48 
and 64 at Mh « 70 GeV 

Alternatively we have used and improved the 
equal weight criterion for the localization of the 
transition. In order to account for mixing of bro- 
ken and symmetric phases at some floating Ph 
inside the metastability range we make the super- 
position ansatz for the histogram: 

Pmip^Pn) = WbPb{p^,PH)+WsPs{p^,pH) + 
WintPint{p^,PH)- (3) 
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For finite volumes there is also a contribution of 
inhomogeneous states containing interfaces. All 
histograms as well as the sum of the appropriate 
thermodynamic weights Wi are normalized. The 
histogram Pm{p'^, Ph) is obtained by reweighting 
towards Ph, merging all data by the Ferrenberg- 
Swendsen method. The pure phase histograms 
should be provided by the same technique from 
data taken outside the metastable range. Our 
data are scarce there, however. 

To obtain the pure phase histograms from runs 
in the metastability region one could use a se- 
lection procedure removing tunneling and excur- 
sions to the "wrong" phase from the Monte Carlo 
history. This method has been used to extract 
pure phase correlation lengths (see below) near T^. 
It is unsuitable for the extraction of the weights 
Wi. We have defined pure phase histograms by 
Gaussian fits to the outer flanks of the full his- 
togram Pm{p^ , Ph), which allows to define the 
weights as functions of Ph and to determine the 
equilibrium point where Ws — wt. More informa- 
tion is given by the probability of mixed phase 
configurations, wtnt = ^ — Wb — Wg at Phc and 
the slopes of Wg and wt at intersection. The infi- 
nite volume limit of the pseudocritical couplings 
is = 0.3435428. We show in fig.|the full his- 
togram Pjn at Pj^^ for the L = 64 lattice and the 
contribution from inhomogeneous (mixed phase) 
configurations. For a lattice of this size (in con- 
trast to the smaller L — 48 or L — 30) the pure 
phase histograms do not overlap, the gap is filled 
exclusively by inhomogeneous states. 



4. Strength of the phase transition 

The present variant of the equal weight method 
for localizing the phase transition has the advan- 
tage to give an immediate estimate for the sur- 
face tension. It is safer than Binder's method 
which uses only the maxima and the minimum 
of Pm{p^, Phc). This is an important improve- 
ment for very asymmetric (typical for Higgs tran- 
sitions) and strongly overlapping histograms (for 
all except the largest lattice sizes). The excess 
free energy per unit area of interphase surfaces a 
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Figure 2. Histograms P„ 
L = 64 at Mh w 70 GeV 



can be obtained from 



and Pint at Phc for 



Wint 
Wb 



exp(- 



(4) 



where the weights are taken at Phc (wb — Ws) and 
A is the surface between the phases. The smallest 
surface A is equal to A = 2L^a^. Thus we obtain 
an estimate, a/T^ « 1.3 x 10"^ at L = 64. 

The latent heat Ae of the transition is given by 
the jump of the average Higgs modulus squared 
according to 



Ae Ml , ^ ^ ^ , 2x 



(5) 



The phase separated histograms with respect to 
at Phc give A(/9^) slightly decreasing with 
growing lattice size. The infinite volume extra- 
polation with gives Ap^ = 0.499(10), which 
amounts to Ae/T^ = 0.024(2). For comparison, 
we have obtained Ae/T^ — 0.20(1) for the case 
Mh = 38.3 GeV (at /3g = 12, too). 
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The equal weight method makes it possible to 
reconstruct the free energy densities of the pure 
phases in the vicinity of the phase equilibrium. 
The latent heat can then be expressed as the jump 
Ae of the energy density by 



Ae 

T4 



Pag' 



dT 

log Ws log Wb 



L3 



L3 



(6) 



with derivatives taken at Phc- Thus we obtain 
Ae/T^ = 0.0238 for L = 64. Both methods give 
compatible results for Ae in the infinite volume 
limit within an l/L^ extrapolation. 

The jump in (p^) at the critical temperature 
(at (3 He) may be translated into continuum units 
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This jump is independent of the 3-dimensional 
renormalization scale /is, and is therefore more 
appropriate to consider than the Higgs conden- 
sate itself. Extrapolating to infinite volume we 
find v/A^ = 0.706 gTc (at /3g = 12). This value 
may be compared with predictions from 2-loop 
perturbation theory (in Landau gauge) Jl3| , pt 
yA^lpert — 0-765 gTc- One has to realize that 
•\/ A(p is not identical with the condensate value 
v{Tc). There is a relation a/A^ = 0.904 u(rc) at 
^3/93 = 0.0957 and with a renormahzation scale 

We conclude that the phase transition at Mh « 
70 GeV shows a clear 2~state signal, and quanti- 
ties like the latent heat, the surface tension and 
the Higgs condensate (all scaled by an appropri- 
ate power of Tc) have a clearly non-vanishing in- 
finite volume limit. For the latent heat and the 
surface tension it must be stressed that, contrary 
to the cases of very small Higgs mass, the transi- 
tion appears to be weaker than predicted by per- 
turbation theory. 



5. The strongly coupled symmetric phase 

This phase is characterized by a mass scale g^T. 
It manifests itself e.g. in the string tension of the 
dimensionally reduced variant of the theory. As 
in the 4-dimensional theory a force can be de- 
fined formally between static non-abelian charges 
through Creutz ratios of Wilson loops. In a pre- 
vious work Q we have calculated the force in the 
case of Mh ~ 35 GeV in the symmetric as well 
as in the broken phase. We have found a string 
tension a — 0.11 (g'^T)'^ in the symmetric phase, 
whereas it vanishes in the broken phase. This 
value has to be compared to cr = 0.13 jg^T)^ for 
the 3-dimensional pure gauge theory [[l5[ . 
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Figure 3. Static force F vs. distance R for Mh 
35 GeV and Mh « 70 GeV 



In fig. ^ we show the force in the symmetric 
phase for the cases Adn ~ 35 GeV and Mh ~ 70 
GeV, both at /3g = 12. From these data for 
Mh ~ 35 GeV the above string tension had to 
be fitted. The errors are large but there is no evi- 
dence for the expected breakdown of confinement 
beyond some screening length. For Mh ~ 70 



6 



GeV (due to the smaller W mass in the symmet- 
ric phase compared to the lighter Higgs case) the 
perturbative contribution extends further in dis- 
tance. Thus we can obtain only an upper bound 
for the string tension, a < 0.1{g^T)^. 

There is another possibility to calculate the 
string tension in the vicinity of the phase transi- 
tion. We consider the ratio of Wilson loops from 
both sides of the phase transition which amounts 
to a subtraction of the perturbative part. This 
allows to anticipate the string tension already at 
much smaller distance. By this method we make 
more precise for Mh ~ 35 GeV the string tension 
a = 0.124 (g2r)2 and predict a = 0.078 (g^T)^ 
for Mh f« 70 GeV. 
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Figure 4. Correlation lengths in Higgs and W 
channel near the transition for Mh ~ 70 GeV 



In order to characterize immediately the sym- 
metric phase we have studied masses in vari- 
ous channels. The behaviour of two correlation 
lengths at Mh ~ 70 GeV with Higgs and gauge 
boson quantum numbers is presented in fig. ^ 
in the vicinity of the phase transition. A discus- 



sion of these results together with the question of 
gauge condensates are postponed to Ref . . 
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